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The critical current of a thin superconducting strip of width W much larger than the Ginzburg- 
Landau coherence length ^ but much smaller than the Pearl length A = 2X^ /d is maximized when 
the strip is straight with defect-free edges. When a perpendicular magnetic field is applied to a 
long straight strip, the critical current initially decreases linearly with H but then decreases more 
slowly with H when vortices or antivortices are forced into the strip. However, in a superconducting 
strip containing sharp 90-degree or 180-degree turns, the zero-field critical current at i? = is 
reduced because vortices or antivortices are preferentially nucleated at the inner corners of the 
turns, where current crowding occurs. Using both analytic London-model calculations and time- 
dependent Ginzburg-Landau simulations, we predict that in such asymmetric strips the resulting 
critical current can be increased by applying a perpendicular magnetic field that induces a current- 
density contribution opposing the applied current density at the inner corners. This efi^ect should 
apply to all turns that bend in the same direction. 

PACS numbers: 74.25.Sv,74.78.-w,74.78.Na 



I. INTRODUCTION 

The critical current of superconductors is of great in- 
terest for both fundamental and practical reasons. On 
the fundamental side, it is of interest to know the precise 
mechanism by which the superconducting state breaks 
down in high applied electrical currents. On the practical 
side, numerous applications of superconductivity require 
the operating current to be as high as possible without 
exceeding the critical current. In this paper we focus 
on thin and narrow strips, where it is known that the 
critical current is dominated by the edge-pinning criti- 
cal current, which is generally much larger than the 
bulk-pinning critical current and can even approach the 
Ginzburg-Landau depairing critical current.^ Pinning in 
thin films by the edge barrier is closely related to pin- 
ning in bulk supercondiictors by the Bean-Livingston 
barrier, ^'^ a barrier against vortex entry that permits 
high currents to be carried by the surface in parallel ap- 
plied magnetic fields. Because in bulk superconductors 
this vortex-entry barrier depends upon both the applied 
field angle and the surface quality, and because there is 
practically no exit barrier, some experiments have found 
that the critical current depends upon the current direc- 
tion, leading to partial rectification'' or diode effect. ^'^ 
Theoretical studies have examined how the quality of the 
surface can alter the Bean-Livingston barrier in bulk su- 
perconductors in parallel fields. Several experimen- 
tal and theoretical studies of thin films^^"^'^ have exam- 
ined the role of edge pinning and found vortex nucleation 
preferentially occurring at defects along the edge or at 



sharp corners where current crowding^*"^° occurs. An 
increase or decrease in the critical current (a diode ef- 
fect) has been found experimentally to depend upon the 
sign of a magnetic field parallel to the surface of a su- 
perconducting strip with a magnetized magnetic strip on 
top, an effect the authors explained chiefly in terms of 
an edge-barrier effect. Roughening one of the edges of 
a straight strip, thereby producing an asymmetric bar- 
rier, has been predicted^^ to lead to rectification upon 
application of a perpendicular magnetic field. 

A recent study^^ has presented calculations of the crit- 
ical currents in thin superconducting strips with sharp 
right-angle turns, 180-degrce turnarounds, and more 
complicated geometries, where all the line widths W were 
much less than the Pearl length^^ A = 2A^/d (film thick- 
ness d, London penetration depth A, < A) but much 
greater than the Ginzburg-Landau coherence length ^. 
That study, in which the critical current was defined as 
the current at which the Gibbs free-energy barrier against 
vortex nucleation is reduced to zero, showed that current 
crowding, occurring whenever the current rounds a turn, 
reduces the critical current below the value it would have 
in a straight strip with defect-free edges. 

In this paper we extend the work of Ref. 23 to consider 
the effect of an applied perpendicular magnetic field H 
upon the critical current. In thin and narrow straight 
superconducting strips for which £^ <^ W <^ A and the 
edge barrier to the nucleation of vortices or antivortices 
dominates the critical current, the critical current Ic{H) 
initially decreases linearly with H When H = 0, the 
critical current is determined by the condition that the 
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current density at the strip's edge reduces to zero the 
Gibbs free-energy barrier against nucleation of a vortex 
on one edge or an antivortox at the opposite edge. When 
H > 0, the current induced by the apphed field is in the 
same direction as the apphed current density at one edge 
but in the opposite direction at the other edge. Thus the 
effect of the apphed field is either to increase the critical 
current for the nucleation of a vortex but to decrease 
the critical current for the nucleation of an antivortex 
or to increase the critical current for the nucleation of 
an antivortex but to decrease the critical current for the 
nucleation of a vortex. In either case, Ic{H) is decreased 
as H increases, because either vortices or antivortices 
are nucleated at a lower value of the applied current. 
The linear decrease with H changes to a slower rate of 
decrease at higher fields, when vortices or antivortices are 
forced to remain in the strip. 

If the superconducting strip is asymmetric, such as 
making a bend to the left, when H = 0, the critical cur- 
rent is determined solely by the condition that the cur- 
rent density around the bend reduces to zero the Gibbs 
free-energy barrier against nucleation of a vortex at the 
inner corner of the bend. When _ff > 0, if the current in- 
duced by the applied field is in the same direction as the 
applied current density at the inner corner, vortex nucle- 
ation occurs at a lower value of the applied current and 
Ic{H) decreases with H. However, if the current induced 
by the applied field opposes the applied current density, 
vortex nucleation occurs at a higher value of the applied 
current and Ic{H) increases with H. 

In Sec. II, we discuss the behavior of the critical sheet 
current Kc{Hz) = Ic{Hz)/W vs applied field in a 
long thin and narrow straight superconducting strip of 
width W . In Sec. Ill, we consider the behavior of long 
thin and narrow strips of width W but with sharp turns 
in the middle. We show that the critical current at zero 
applied field is reduced below the corresponding critical 
current of a straight strip with no turn, but we predict 
that its critical current can be increased by applying a 
magnetic field of the right polarity. In Sec. IV, we re- 
port simulations using time-dependent Ginzburg-Landau 
equations that generally confirm these predictions. In 
Sec. V, we summarize our results and discuss possible 
applications in superconducting nanowire single-photon 
detectors. The Appendix contains details of how to cal- 
culate the functions needed in Sec. III. 



II. Ic{H^) FOR A STRAIGHT STRIP 

Consider a superconducting strip of London penetra- 
tion depth A, thickness d {d < A), two-dimensional 
screening length (Pearl length'^^) A = 2X^/d, width 
W {W ^ A), and Ginzburg-Landau coherence length 
<C W), centered on the xy plane in the region 
< y < W. A current I flows in the x direction. Be- 
cause <C A, the corresponding sheet-current density 
Kj = xKj = xl /W is very nearly independent of y?^ In 



the presence of an applied magnetic field H = zH^, the 
magnetic field penetrates very nearly uniformly through 
the strip. The corresponding iJ-induced sheet-current 
density Kh = xKhx (y) can be calculated from the Lon- 
don equation K = —{2/fioA)[A + (0o/27r)V7], where A 
is the vector potential {B = V x A), <j)o = h/2e is the 
superconducting flux quantum, and 7 is the phase of the 
order parameter: KHx{y) = (2i?^/A)(y - T4^/2). 

We next use a London-model description of vortices to 
estimate the critical current when the critical current due 
to the edge barrier greatly exceeds the critical current due 
to bulk pinning. We start by writing down the Gibbs free 
energy G of a vortex nucleating at y = T4^ or an antivor- 
tex nucleating at y = 0, denoting 5 as the distance from 
the edge {5 = W — y iov the vortex, which carries mag- 
netic flux in the z direction, and 5 = y for the antivortex, 
which carries magnetic flux in the —z direction) 

4>l , \2W . f7rS\i 

T'^SiW-S), (1) 

where the flrst term is the self-energy of the vortex or 
antivortex, accounting for its interactions with an infinite 
set of images, the second term is the negative of the work 
done by the source of the applied current as the vortex or 
antivortex moves from the edge to the coordinate 5, and 
the third term is the negative of the work done by the 
source of the applied magnetic field as the vortex (upper 
sign) or antivortex (lower sign) moves from the edge to 
the coordinate 5. Note that G is minimized for the vortex 
when Hz is positive and for the antivortex when is 
negative. 

As in Ref. 23, to estimate the Gibbs free energy barrier 

in different geometries, in Eq. (1) and later in this paper 
we use the simplifying assumptions of a London-model 
vortex and its neglect of the vortex-core energy, which 
are known to lose accuracy when the vortex is close to 
the sample edge. Although the numerical factors in our 
expressions for the critical current therefore are probably 
not accurate, we expect the qualitative behavior of their 
geometry dependence to be correct. 

A. Linear behavior for small H 

To calculate the nucleation of a vortex or antivortex 
when ^ <C W , we need to examine the behavior of G only 
for small values of 5, for which, to good approximation, 

G = In (^) - cl^oKid T MHzW/K)5, (2) 

where the upper (lower) sign holds for the nucleation of 
a vortex (antivortex). The free-energy barrier occurs at 
5 = Sb- Setting dG/dS = there, we obtain 



2niJ,oA{Ki ± HzW/A) ' 
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FIG. 1: Critical sheet current for vortices KcviH^) and an- 
tivortices Kca{Ha) (normalized to Kca) vs Hz (normalized to 

J^stop). 



which describes a force balance between the repulsive 
Lorentz force (poiKi ± H^W/A) and the attractive force 

of the nearest image (/)Q/27r/ioAi5. Setting G = a.t Sf, 



yields 5b = Sc = 
sheet current Kj 



= 



1.36^ <C W and the critical 
where 



(4) 
(5) 



e = 2.718... is Eulcr's number and Kcs denotes the criti- 
cal sheet current for a long straight strip in the absence of 
an applied field. Note that for the nucleation of a vortex 
(antivortcx) the critical sheet current is decreased (in- 
creased) for positive Hz , because at the nucleation point 
the iJ-induced current density is in the same (opposite) 
direction as the applied current density. 

As the applied sheet current Kj is increased, a volt- 
age appears along the length when the smaller of the 
two critical sheet currents in Eq. (4) is reached. For 
Hz > {Hz < 0) the critical current is reached when 
conditions are favorable for vortices (antivortices) to be 
nucleated. Because of this symmetry, with identical bar- 
riers for vortices or antivortices on opposite sides of the 
strip, the critical sheet current for small H = \Hz\ is 
therefore Kc{H) = Kcs ~ HW/A. This result also can 
be written as Kc{H) = Kcs[l - {eW/2i){H/Hc2)], where 
H^2 = <^o/27rMo€'- See Fig. 1. 



B. Behavior for large H 

The linear decrease of Kcv{Hz) = - H^W/A for 
the nucleation of vortices (hence the subscript v) given 
in Eq. (4) applies only for relatively small values of Hz- 
At the critical current the nucleating vortex must be able 
to travel all the way across the strip and annihilate with 
its image on the opposite side. However, it can be seen 
from Eq. (1) that large values of Hz can produce a free- 
energy minimum for 6 < W . Accordingly, as Hz in- 



creases, there is a special value of Hz , which we denote as 
Hstop, at which this free-energy minimum first appears. 
A nucleating vortex therefore comes to a stop at the dis- 
tance (^niin where the minimum occurs. When £,/W <C 1, 
W{1 - ^e£,/2W), which is very close to the edge 
where the vortex would have annihilated, and the value 
of ilstop is given to good approximation by 



-^stop — 2W 



H 



c2- 



(6) 



The linear decrease of Kcv{Hz) therefore ceases at Hz 
-ffstop) where (to good approximation when ^ ^ W) 

Kcv{Hstop) = Kcs/2. 



(7) 



For Hz > Hstop, Kcv decreases more slowly with Hz, 
because the vortices that have stopped within the strip 
generate an additional current density opposing the ap- 
plied current density where vortices are nucleated. To 
calculate Kcv (Hz) accurately becomes a complicated nu- 
merical problem, because one must take into account the 
effect of all the vortices that temporarily reside in their 
local Gibbs free-energy minima. 

A crude estimate of the field dependence of Kcv{Hz) 
for Hz > Hstop can be obtained using an approach analo- 
gous to that used in Ref. 29, in which the field-dependent 
critical current due to surface barriers was calculated for 
superconducting strips of width W ^ A, the limit oppo- 
site to that of interest to us here. We consider the behav- 
ior just below the critical current in a field Hz > i?stop 
when an array of vortices is held in place by the applied 
field. Our approximation is to replace the discrete vortex 
array by a stationary distribution of vortices of uniform 
density Uy = ^^Hz/ipo within a band < y < y^, where 
yy = W{1 — Hstop/Hz). From the London equation and 
the fact that a vortex contributes no net current along 
the sample length, it follows that this band of vortices 
generates a current density 



Ky,{y) = {Hz/WA)[2W{yy-y)-ytl < y < (8) 
= -HzyllWA, yv<y<W.{Q) 

The net sheet-current density in general is 

= Ki + KH^{y) + Ky^iy), (10) 

but since we must also have K^ {y) = within the vortex- 
filled region {0 < y < y^) so that these vortices do not 
move, we obtain 



Ki = {Hz/WA){W -y2f. 



(11) 



At the critical current, the vortices present in the 
region < y < y-z produce the additional sheet cur- 
rent Kyx{W) = —Hzyl/WA. To determine the condi- 
tion for which a new vortex is nucleated using Eq. (2), 
we must therefore add to the right-hand side an addi- 
tional term (paHzy^S/WA. The quantity in parenthe- 
ses in the denominator of Eq. (3) is then replaced by 
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Ki + H^W/K - H^yl/WK. When this quantity reaches 
the value K^s, Kj reaches its critical value, Kcv{Hz), 
such that 

Ki = - H,W/A + H,yl/Wk. (12) 
Setting Ki = Kcv{Hz) in Eqs. (11) and (12), we obtain 

for H:, > -ffstop, 

2/2 = W{l-H,top/H,), (13) 
K,^iH,) = KM{Hsiop/2H,). (14) 

Thus, in this approximation, the linear behavior 
K,,{H,) = K,{0){1 - HJ2H,,,^) for < J?, < F.top 
turns into a slower Kcv{Hz) = Kc{0){Hstop/'^Hz) de- 
crease for Hz > -ffstop with no change in slope at = 
iJstop, a behavior similar to that found in Ref. 29. 

Note that the calculations given here assume that the 
sheet-current density in the strip is much higher than 
the critical sheet current attributable to bulk pinning of 
vortices or antivortices. The effects of bulk pinning could 
be accounted for using a treatment similar to that in Ref. 
3. 

The solid curves in Fig. 1 show the critical sheet cur- 
rent of a long straight strip (normalized to Kcs) vs Hz- 
For Hz > 0, the critical sheet current is KcviHz), which 
is determined by the nucleation of vortices (hence the 
subscript v) at y = W, while for Hz < 0, the critical 
sheet current is Kca{Hz), which is determined by the nu- 
cleation of antivortices (hence the subscript a) at y = 0. 
The dashed curves show the extension of the linear por- 
tions of the curves. The mirror symmetry about the ver- 
tical axis at Hz occurs because of the mirror symmetry 
of the sample geometry about y — W/2. As we show in 
the next section, if the sample geometry is asymmetric, 
there is a dramatic difference in the behavior of the crit- 
ical sheet current depending upon the direction of the 
applied field. 



III. /c(Hz) FOR STRIPS WITH LEFT TURNS 

Next wc consider strips of width W carrying cur- 
rent KjW making left turns around a sharp 90-degree 
left turn [Fig. 2(a)], a sharp rectangular 180-degree 
turnaround [Fig. 3(a)], or a sharp 180-dcgrcc turnaround 
[Fig. 4(a)]. In Ref. 23, using a London-model description 
of vortices and antivortices, Clem and Berggren devel- 
oped a procedure for estimating how the critical sheet- 
current density in the absence of an applied field Kc{0) is 
reduced because of current crowding at the inner corners 
of turns and bends. For each sample geometry consid- 
ered, they carried out a calculation of the Gibbs free en- 
ergy similar to that in Eq. (2) and showed that when 
^ <^ W, the critical sheet-current density can be ex- 
pressed as Kc{0) = KcsR, where Kcs is the critical sheet 
current for a long straight strip in the absence of an ap- 
plied field [Eq. (5)] and i? is a critical-current reduction 
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FIG. 2; (a) Current flow in a strip carrying current KiW 
around a left-hand 90° turn, shown by the contour plot of 
its stream function Siz{x,y). The contours correspond to 
streamlines of the induced sheet-current density Ki{x, y), and 
the arrows show the direction of the current. Note the current 
crowding that occurs near the inner corner E. (b) Clockwise 
current flow induced by a positive applied flold Hz, shown 
by the contour plot of the stream function SHz{x,y). The 
contours correspond to streamUnes of the fleld-induced sheet- 
current density KH{x,y), and the arrows show the direction 
of the current. Note the current crowding that occurs near 
the inner corner E. 



factor (0 < i? < 1). For left-hand turns, we can ap- 
pend the subscript v to the critical sheet-current density 
Kcvi^) as a reminder that the critical current occurs at 
the threshold for the nucleation of vortices at the sharp 
inner corners. Antivortices could be nucleated along the 
long, straight portions of the outer boundaries only at a 
higher current density, ifca(O) = i^cs- 

In this section, we extend the approach of Ref. 23 to 
study the effect of an applied perpendicular magnetic 
field H = zHz upon the critical sheet-current density 
for the sample geometries shown in Figs. 2-4. When 
the applied current makes only left turns, as shown in 
Figs. 2(a), 3(a), and 4(a), the critical current for small 
Hz is dominated by the nucleation of vortices at the 
sharp inner corners, and we therefore seek expressions for 
Kcv{Hz), where the subscript v denotes vortices, which 
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FIG. 3: (a) Current flow in a strip carrying current KjW 
around a sharp rectangular 180° turn, shown by the contour 
plot of its stream function Siz{x, y). The contours correspond 
to streamlines of the induced sheet-current density Ki{x,y), 
and the arrows show the direction of the current. Note the 
current crowding that occurs near the inner corners E and 
F. (b) Clockwise current flow induced by a positive applied 
field Hz, shown by the contour plot of the stream function 
SHz{x,y). The contours correspond to streamlines of the 
field-induced sheet-current density KH{x,y), and the arrows 
show the direction of the current. Note the current crowding 
that occurs near the inner corners. 



carry magnetic flux in the +z direction. The critical- 
current calculation requires knowledge of the detailed 
form (especially the behavior of the divergences near 
the inside corners) of the applied sheet-current density 
Ki{x,y) = V X Si{x,y), where Si{x,y) = zSi;,{x,y) 
is its stream function, and the magnetic-field-induced 
sheet-current density Kh = V x Sh, where Sh = zShz 
is its stream function. Figures 2(a), 3(a), and 4(a) show 
contour plots of the stream functions Sj{x,y), and Figs. 
2(b), 3(b), and 4(b) show contour plots of the stream 
functions Sh{x, y). For each of the three geometries, the 
field-dependent critical sheet-current densities for the nu- 
cleation of vortices at the inner corners can be expressed 
as 

K,,{H,) = K,,{0) - aWH,/A, (15) 
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FIG. 4: (a) Current flow in a strip carrying current KjW 
around a sharp 180° turn, shown by the contour plot of 
its stream function Siz{x,y). The contours correspond to 
streamlines of the induced sheet-current density Ki{x, y), and 
the arrows show the direction of the current. Note the cur- 
rent crowding that occurs near the origin O. (b) Clockwise 
current flow induced by a positive applied held Hz, shown 
by the contour plot of the stream function SHz{x,y). The 
contours correspond to streamlines of the fleld-induced sheet- 
current density KH{x,y), and the arrows show the direction 
of the current. Note the current crowding that occurs near 
the origin. 



where Kcv{0) = KcsR- The specific forms of the critical- 
current-reduction factor R and the field-slope parameter 
fj (0 < (7 < 1) depend upon the sample geometry. The 
details of the calculations of Si{x, y), Sh{x, y), R, cr, and 
Kcv{Hz) are given in the Appendix. 

Note that Kcv decreases with increasing H^, because 
the field-induced sheet-current density is in the 

same direction as the applied sheet-current density at 
the vortcx-nucleation site. Kcv increases for negative Hz 
of increasing magnitude, because the field-induced sheet- 
current density is then in the opposite direction as 
the applied sheet-current density at the vortex-nucleation 
site. 

We next examine how the application of a magnetic 
field H = zHz affects the critical current due to the 
nucleation of antivortices at the outer boundaries of 
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FIG. 5: (a) Left turns: Generic behavior of the critical sheet 
current for vortices Kcv(Hz) and antivortices Kaa{Ha) (nor- 
malized to Kcs) vs Hz (normalized to A.Kcs/W), illustrated 
here for a sharp 90-degree turn when the current / flows 
down the vertical portion of the strip shown in Fig. 2(a) 
and makes a left turn axound the corner E. (b) Right turns: 
Same as (a) except for the case that the current / flows in 
through the horizontal portion of the strip shown in Fig. 
2(a) and makes a right turn around the corner E. The val- 
ues £,/W = 0.047 and R = 0.5 were assumed for both plots, 
such that WHp/h.Kcs = 0.287 and K^p/Kcs = 0.713. 



the structures shown in Figs. 2-4. This occurs where 
the magnitude of the sheet-current density K{x,y) = 
Ki{x,y) + Knix^y) has its maximum value. It can be 
shown that, except for negative values of of large mag- 
nitude, this maximum occurs on the outer boundary at 
large distances from the outer corners, where the spatial 
dependence of K{x,y) = Ki{x,y) + KH{x,y) becomes 
the same as in a long, straight strip of width W. Thus, 
the critical sheet-current density at which antivortex nu- 
cleation occurs is 



(16) 



The critical current for the nucleation of antivortices in- 
creases linearly from the value Kcs for > G and de- 
creases linearly for < 0, just as in the case of an 
infinitely long strip, as given in Eq. (4) and shown in 
Fig. 1 by the curve Kca{Hz)- 

The solid lines in Fig. 5(a) show the generic behavior 
of the critical sheet-current density for all three left-hand 
turns shown in Figs. 2-4, but illustrated using parameters 
for the sharp 90-degree turn, R = 0.5 and a = SG/tt^ = 
0.742, where G = 0.915965... is Catalan's constant. The 



overall critical current of films with sharp left-hand turns 
is determined by the critical current for the nucleation of 
vortices for all values of > H^ax a-nd by the critical 
current for the nucleation of antivortices for all values of 
Hz < Hjnax, where Hmax is the applied field at which the 
critical currents for vortices and antivortices are equal. 
From Eqs. (15) and (16) we see that when the applied 
field is Hmax = —Hp, where 



AKcs\{l-R) 2^(1 -i?) 



(l + a) eW{l + a] 



H, 



c2, 



(17) 



the maximum overall critical sheet current has the peak 
value Kc{Hmax) = Kcp, where 



Kcv{0) R{l + ay 



(18) 



The latter equality in Eq. (17) arises from Eq. (5) and the 
well-known expression iiqHc2 = ^o/27r^^. In summary, 
we see that by applying a negative applied field H^ = 
Hmax = ^Hp, the critical current, which in zero field is 
suppressed by current crowding at the sharp inner corner, 
can be increased by the factor (i? -|- a)/R{l + a). This 
factor is very large when R is small, is equal to 1.43 for 
the sharp 90-degree turn when R = 0.5 and a = 0.742, as 
shown in Fig. 5, and approaches 1 as — >■ 1. Expressions 
for i? and a for all three geometries shown in Figs. 2-4 
are given in the Appendix. 

The generic behavior of the field-dependent critical 
current when the current makes right turns is shown in 
Fig. 5(b). Consider what happens in a strip with a sharp 
90-dcgree turn as in Fig. 2(a) when the direction of the 
current / is reversed. At if^ = 0, the critical current is 
now determined by the onset of antivortex nucleation at 
E. This critical current increases linearly with H^, but 
the maximum overall critical current is reached when 
Hmax = Hp and the critical current for antivortex nu- 
cleation at E becomes equal to the critical current for 
vortex nucleation along one of the straight sides far from 
the outer corner B. The overall critical sheet current is 
shown by the solid lines in Fig. 5(b). 

The above analysis, which assumes ^ <C W, is in- 
tended to describe the critical current only at relatively 
low fields, where the behavior is linear in the applied field. 
At high applied fields, vortices or antivortices are forced 
into the strip, and the dependence of the critical current 
becomes nonlinear, similar to the behavior discussed in 
Sec. II B. The effect of the vortex distribution inside the 
sample on the magnetic-field dependence of the critical 
current will be addressed in the next section. 



IV. NUMERICAL SIMULATIONS WITHIN THE 
TIME-DEPENDENT GINZBURG-LANDAU 
THEORY 

In order to confirm the theoretical predictions in the 
preceding sections, we performed numerical simulations 
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within the time-dependent Ginzburg-Landau (GL) the- 
ory. We consider superconducting strips (with thick- 
ness d X and width w ^ A = 2X/(P) with muhiple 
sharp turns in the presence of a transport current (ap- 
pUed through the normal contacts) and a perpendicular 
magnetic field of magnitude H (see Fig. 6). 



and chose A = iiQ{—Hy/2,Hx/2). We solve the above 
coupled non-linear differential equations self-consistently 
in 2D using Euler [for Eq. (19)] and multi-grid^^ (for the 
equation of the electrostatic potential) iterative proce- 
dures. 



A. Theoretical approach 



Superconducting strip with a 90-degree turn 



To study the dynamics of the superconducting con- 
densate in such complex structures, we used the follow- 
ing generalized time-dependent Ginzburg-Landau (GL) 
equation' 



30. 



d_ 

di 



-lif- 



t/j^iV - iA) V 

+(i-i^nv^,(i9) 



which is coupled with the equation for the electrostatic 
potential Aip = div{Im[-0*(V — iA)-0]}. Here distance is 
scaled to ^, the vector potential A is in units of (/)o/27r^, 
time is in units of the GL relaxation time Igl — Mo-^^/Pn 
{pn is the normal-state resistivity), and voltage is scaled 
to Vq = 0o/27rtGL- The coefficient u, which governs the 
relaxation of the order parameter (i.e., the ratio between 
relaxation times for the phase and the amplitude of ?/;), 
and the material parameter 7 are chosen as u = 5.79 
and 7 = 10, which are found within the microscopic BCS 
theory for superconductors with weak depairing.'^'^ Us- 
ing the normal-state resistivity p„ = 18.7 /iilcm, zero- 
temperature coherence length ^(0) = 10 nm and pene- 
tration depth A(0) = 200 nm, which are typical for Nb 
thin films, "^^ one can obtain toL ~ 2.69 ps and Vq ~ 0.12 
mV near Tc- We use superconducting- vacuum boundary 
conditions (V — iA)tp\n = and V(/3|„ = at all sam- 
ple boundaries, except at the current contacts where we 
use ■0 = and Vipln = ~j, with j being the applied 
current density in units of jo = <^o/27r/zo^A^. Assuming 
that W <^ A = 2X/(P we neglect demagnetization effects 




FIG. 6: (Color online) The model systems: Superconducting 
strips (of length L, width W and thickness d) with (a) 90° and 
(b) 180° turns. The current is applied through normal-metal 
contacts (and flows always such that it makes a left turn) and 
the output voltage is measured at a small distance away from 
these leads. The magnetic field is applied perpendicular to 
the sample either in the z- (positive field) or the -z-direction 
(negative field). 



We begin our analysis by demonstrating the properties 
of a superconducting strip with a 90° turn [see Fig. 6(a)] 
by constructing the time-averaged voltage vs applied cur- 
rent (/ — V) and the voltage vs time V{t) characteristics 
of the sample for a constant magnetic field. The time- 



0.10 




FIG. 7: (Color online) Current-voltage characteristics of the 
sample with a 90° turn for different values of (a) negative 
and (b) positive magnetic field. Right axis in (a) shows the 
voltage values for — — 0.025-ffc2- Dimensions of the sample 
are L — 32^ and W = 15.5^. Insets show snapshots of the 
Cooper-pair density at the current and magnetic field values 
indicated by open circles in the I — V curves. Panels 1-6 show 
snapshots of the Cooper-pair density (1, 4), the phase of the 
order parameter (2, 5), and the supercurrent distribution (3, 
6) for vortex nucleation at zero magnetic field (1-3) [point 
"c" in (a)] and antivortex nucleation at ~ ~0.05Hc2 (4- 
6) [point "d" in (a)]. White arrows indicate the direction of 
vortex or antivortex motion. 
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averaged voltage shown in Figs. 7, 11, and 13 is equiv- 
alent to V = (pQV, where u is the net rate with which 
vortices cross the strip in one direction and antivortices 
cross the strip in the opposite direction, in accord with 
the Josephson relation. '^'^'■^^ As a representative example, 
we consider a superconducting strip with length L = 32^ 
and width W = 15.5^, the I —V characteristics of which 
are shown in Fig. 7 for different values of negative (a) and 
positive (b) magnetic field. We first discuss the results 
for negative direction of the magnetic field, which is the 
situation when the screening (Meissner) currents oppose 
the applied current near the sharp inner corner of the 
sample. In the absence of the magnetic field [solid black 
curve in Fig. 7(a)] zero resistance of the sample is main- 
tained up to a threshold current density = 0.3275jo, 
above which the system goes into the resistive state with 
a finite- voltage jump. This resistive state is characterized 
by the periodic nucleation of vortices near the inner cor- 
ner where the current density is highest (see panels 1-3 
in Fig. 7). This vortex is driven further by the Lorentz 
force towards the outer corner of the sample where it 
leaves the sample (see panel 1). The nucleation rate of 
vortices in the inner corner increases with further increas- 
ing the applied current, and at sufficiently large currents 
the system transits to a higher dissipative state, char- 
acterized by fast-moving (kinematic) vortices. '^^■'^^ The 
critical current of the sample jc considerably increases 
with applying small negative magnetic field [dashed red 
curve in Fig. 7(a)]. This is because the Meissner currents 
reduce the current crowding at the sharp inner corner. 
However, at larger values of the negative field the critical 
current becomes smaller [see dot-dashed green curves in 
Fig. 7(a)], because antivortices start penetrating the sam- 
ple (see panels 4-6 in Fig. 7). These antivortices (com- 
pare phase plots in panels 2 and 5) nucleate away from 
the corners of the sample (where the current density is 
maximal) and are driven across the sample (white arrows 
indicate the direction of motion). Thus, the critical cur- 
rent of the sample has a non-monotonic dependence on 
the negative magnetic field. At intermediate values of 
the magnetic field, we observed the coexistence of vor- 
tices and (kinematic) antivortices as shown in the inset 
of Fig. 7(a) (see also discussion of Fig. 10). 

The dashed red curve in Fig. 7(b) shows the I — V 
curve of the sample for positive magnetic field. For this 
direction of the field the Meissner currents add to the 
applied current near the inner corner of the sample, re- 
ducing the surface barrier for the nucleation of a vortex 
there. Indeed, analysis of the temporal characteristics of 
the sample (not shown here) shows that the vortices al- 
ways nucleate near the inner sharp corners, as illustrated 
in the inset of Fig. 7(b). Thus, the critical current is 
solely determined by the vortex entry at the inner cor- 
ner, and it becomes a monotonically decreasing function 
of the magnetic field. 

Figure 8 summarizes our findings, where we plot the 
resistive-state transition current jc as a function of ap- 
plied magnetic field (see top axis for the fiux going 
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FIG. 8: (Color online) The critical current jc of the sample 
in Fig. 7 as a function of negative (filled black circles) and 
positive (open red circles) applied magnetic field. The results 
are normalized to the critical current at zero magnetic field. 
Top axis shows the flux (in units of the flux quantum (/)o) 
through the sample. Panels 1-4 show contour plots of for 
magnetic field values indicated on the jc{H) curves and for 
current values just below jc- 

through the sample). At relatively low negative fields 
(filled black circles), the critical current increases with 
increasing magnetic field; up to 30% enhancement can be 
achieved for the given parameters of the sample. How- 
ever for high magnetic fields, jc decreases again with 
increasing negative field, because magnetic-field-induced 
antivortices start to nucleate [see the discussion of Fig. 
9(b)]. At larger fields, the critical current shows a 
diffraction-like pattern as a function of magnetic field. 
This is because different antivortex states are stabilized 
before the system transits to the resistive state (see pan- 
els 1 and 2). For the positive direction of the magnetic 
field (open red circles), the critical current is a monoton- 
ically decreasing function of Hz', at relatively low fields 
jc has a linear dependence on the field. At higher fields, 
vortices start penetrating the sample (see panels 3 and 
4) and the jc{H) curve becomes nonlinear - the behavior 
found for straight superconducting strips (see the discus- 
sion in Sec. IIB and also Ref. 3). Notice that the critical 
current for negative magnetic field is always larger than 
the one for positive field due to the suppression of current 
crowding at the sharp inner corner. 

To get a better insight into the dynamics of vortices 
in the system, we plotted in Fig. 9(a) the time evolution 
of the output voltage at zero magnetic field and for an 
applied current just above jc, together with snapshots 
of the Cooper-pair density at times indicated in the V{t) 
curves. The output voltage oscillates periodically in time 
(one full period is shown) with a minimum correspond- 
ing to the Meissner state [inset 1 in Fig. 9(a)]. This 
voltage can be described as a periodic sequence of single- 
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FIG. 9: (Color online) Voltage vs time response of the sample 
in Fig. 7 at j = 0.34jo for (a) ^ and (b) = -0.05ii'c2, 
upon transition to the resistive state [see points "c" and "d" in 
Fig. 7(a)]. Insets show snapshots of the Cooper-pair density 
at the times indicated on the V{t) curves, which illustrate the 
motion of vortices (a) and antivortices (b). 



vortex pulses, each of integrated area J Vdt = (/iq.^^'^® 
With time, a vortex nucleates at the inner corner of the 
sample (inset 2) where the current density is higher due 
to the current crowding. This nucleation process cor- 
responds to a maximum in the V(t) curve. By means 
of the Lorentz force this vortex moves deeper inside the 
sample towards the outer corner of the sample (inset 3). 
Note that the voltage decreases during the motion of the 
vortex, indicating that the vortex slows down during its 
motion. The expulsion of the vortex also leads to an ex- 
tra voltage peak (inset 4). Note that for this value of the 
current there is only one vortex in the sample at a given 
time; however, for larger currents more than one vortex 
can be present [see the inset in Fig. 7(b)]. 

Figure 9(b) shows the voltage characteristics of the 
sample together with the evolution of the vortex state 
for negative magnetic field at = — 0.05i?c2- The 
output voltage shows periodic oscillations with several 
maxima and minima. The global minimum corresponds 
to the Meissner state (point 1 and inset 1). At a later 
time two antivortices penetrate the sample, one after the 
other (insets 2 and 3), leading to local maxima in the 
V{t) curve just to the left of points 2 and 3. Note that 
antivortex nucleation does not occur at the outer cor- 
ner of the sample because the current is minimal there. 
Rather they nucleate away from this corner, where the 
current density is maximal (see panel 6 in Fig. 7 for the 



current distribution in the sample). These antivortices 
temporarily reside inside the sample, which leads to lo- 
cal minima in the V{t) curve to the right of points 2 
and 3. As time goes on, the antivortices leave the sam- 
ple, one after the other (insets 4 and 6), producing local 
maxima in the V(t) curve (points 4 and 6), beyond which 
the superconducting condensate relaxes towards its ini- 
tial state. After some time, a new antivortex penetrates 
the sample and the entire antivortex entry-exit sequence 
repeats. The time-dependent voltage shown in Fig. 9(b) 
thus can be regarded as a periodic sequence of double 
pulses produced by two antivortices crossing at nearly 
the same time. The integrated area of each double pulse 
is / Vdt — 2cf)Q. Overall, for larger values of the negative 
field the resistive state is characterized by the periodic 
entrance of antivortices. 

Figure 10 (a) shows the V{t) characteristics of the sam- 
ple for Hz = — 0.025-ffc2, where the critical current of the 
sample reaches its maximal value (i.e., jc = 0.425jo). 
Surprisingly, the resistive state in this case is charac- 
terized by the motion of both vortices and antivortices, 
which move in opposite direction without annihilation 
(see white arrows in panel 1 in Fig. 10 for their direction 
of motion). Vortices nucleate at the inner corner of the 
sample periodically in time and move along the diago- 
nal direction (panels 2-5), whereas antivortices nucleate 
at the outer edge of the sample and move across it, as 
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FIG. 10: (Color online) (a) Time evolution of the voltage of 
the sample in Fig. 7 for j = 0.425jo and = — 0.025//c2- 
Panels 1-5 show the time evolution of the superconducting 
condensate at the time intervals indicated on the V{t) curve, 
(b) The Fourier power spectrum of V{t). The nucleation rates 
of vortices and antivortices are indicated. 
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FIG. 11: (Color online) I — V characteristics of the sample 
(L = 32^ and W = 15.50 with a 180° turnaround for different 
values of (a) negative and (b) positive magnetic field. Panels 
1-4 show snapshots of li/ip at field and current values indicated 
on the I — V curves. White arrows indicate the direction of 
vortex motion. 



we discussed previously in connection with Fig. 9. How- 
ever, antivortices move much faster than vortices, they do 
not retain their circular shape, and they create a channel 
with suppressed order parameter. This is due to the more 
uniform current flow in the regions away from the corner, 
which was shown to be the condition for the formation 
of fast-moving kinematic vortices. The speed of these 
kinematic vortices decreases with further increasing the 
applied magnetic field. Coexistence of fast- and slow- 
moving vortices has also been found in straight supercon- 
ducting samples. Thus, the output voltage is character- 
ized by fast oscillations due to the kinematic vortices, the 
amplitude of which is periodically modulated due to the 
slow motion of vortices, as seen from the Fourier power 
spectrum of the V(t) curve [see Fig. 10(b)] 



Superconducting system with a 180-degree 
turnaround 



In what follows, we study how the critical current 
of a superconducting strip with a sharp 180-degree 
turnaround [see Fig. 6(b)] is affected by the external field. 
As our main results, we plotted in Fig. 11 the calcu- 
lated I—V characteristics of the sample with dimensions 
W = 15.5^ and L = 32^ (i.e., the gap separating two 
parts of the sample is A = 1^) for different values of neg- 
ative (a) and positive (b) values of the applied magnetic 
field. As in the case of the sample with a 90-degree turn. 



vortices are bound to nucleate at the sharp inner corner in 
the absence of a magnetic field (see the solid black curve 
and panel 1 in Fig. 11), the motion of which determines 
the resistive-state transition current. This critical cur- 
rent increases with increasing negative field until some 
threshold field (see dotted blue curve), above which jc 
decreases again (dot-dashed green curve). At these val- 
ues of the field, antivortices penetrate the sample away 
from the corners (panel 3), leading to energy dissipation 
in the system. The coexistence of vortices and antivor- 
tices is also observed at intermediate values of the field, 
as shown in panel 2 of Fig. 11. For positive direction of 
the magnetic field, vortices preferentially nucleate at the 
inner corner (see panel 4) where the screening currents 
add to the applied current. Thus, the critical current 
becomes a decreasing function of the field, because vor- 
tex entries are shifted to lower currents as compared to 
the zero-applied-field case. Note that the critical current 
density jc is always smaller than the one obtained for 
the 90-degree turn sample due to larger current crowd- 
ing near the inner corners (compare the I ^ V curves in 
Figs. 7 and 11). 

To determine how the critical current is affected by 
the applied field and, in particular, to show to what ex- 
tent the critical current of the sample can be increased 
by the negative magnetic field, we present in Fig. 12 the 
critical current of the system jc as a function of the mag- 
netic field. jc{H) curves show behavior similar to that 
in the case of a single 90-degree turn (see Fig. 8): (i) 
jc increases with increasing negative field (filled black 
circles) until some threshold field, beyond which jc de- 
creases again due to the penetration of antivortices; (ii) 
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FIG. 12: (Color online) The critical current jc/jc{0) of the 
sample in Fig. 11 as a function of negative (filled black cir- 
cles) and positive (open red circles) applied magnetic field. 
Panels 1-4 show contour plots of \'ip\'^ for magnetic field val- 
ues indicated on the jc{H) curves and for current values just 
below jc- 
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FIG. 13: (Color online) (a) I — V characteristics of the sample 
with L = 32^ and W = 13.5^ (i.e., S = 5^ with a 180° 
turnaround for different values of the negative field. Insets 
show snapshots of \ip\'^ at field and current values indicated 
on the I — V curves, (b) jc{H) curves of the sample for 
negative (solid black circles) and positive (open red circles) 
magnetic field. 



for positive direction of the field (open red circles), jc lin- 
early decreases with field because of the increased current 
density at the corners, which reduces the energy barrier 
for the nucleation of vortices; (iii) at larger field values 
different vortex (antivortex) patterns are stabilized in the 
system (see panels 1-4 in Fig. 12) which results in a non- 
linear dependence of the critical current on both negative 
and positive magnetic fields. 

Finally, we discuss the effect of the opening A of 
the sample on our findings [see Fig. 6(b)]. Figure 13 
(a) shows the I ~ V characteristics of the sample with 
L = 32^ and A = 5^ for several values of the negative 
magnetic field. At zero magnetic field (as well as for 
positive applied fields) vortices nucleate at the two in- 
ner corners of the sample, one after another [see inset 1 
in Fig. 13(a)], and move along the diagonal direction, as 
in the case of the 90-degree turn sample. The zero-field 
critical current of the sample is considerably larger than 
the one for the sample with a narrow gap (see Fig. 11), 
because of the smaller current crowding near the inner 
corners, jc increases with increasing negative magnetic 
field [filled black circles in Fig. 13(b)] until some thresh- 
old value, above which jc is reduced due to the formation 
of antivortices [inset 3 in Fig. 13(a)]. As before, jc is a 
decreasing function of the positive magnetic field. 

In conclusion, our numerical simulations confirm the 
following theoretical findings for superconducting sam- 
ples with sharp (90- or 180-degree) turns: (i) in the ab- 
sence of the magnetic field, vortices nucleating at the 
inner sharp corners determine the resistive-state transi- 



tion current, which is smaller than the one for a straight 
strip with no turn due to the current crowding at the 
inner corners; (ii) the critical current can be increased 
by applying a perpendicular magnetic field of the right 
polarity, because of the reduction of the applied current 
density at the inner corners by the screening currents; 
(iii) for the other direction of the magnetic field, the crit- 
ical current decreases linearly with field, because at the 
nucleation point the field-induced current is in the same 
direction as the applied current; (iv) for larger values of 
the magnetic field the critical current decreases (but with 
slower rate) with the field in either direction; (v) oscilla- 
tions in the jc{H) curves are observed at larger fields as a 
consequence of different vortex states, before the system 
transits to the resistive state. 



V. SUMMARY AND DISCUSSION 

In this paper we have considered thin (thickness d < X) 
and narrow (width W <^ A = 2\^/d and £,<^W) super- 
conducting strips with sharp turns in the middle. Using 
a London-model description of vortices, we showed the- 
oretically that when an applied current / flows through 
the strips shown in Figs. 2-4 and makes left turns, the 
critical sheet current, dominated by the onset of vortex 
nucleation at the sharp inner corners, is suppressed by 
current crowding. We have shown that this critical cur- 
rent further decreases linearly as a positive magnetic field 
Hz is applied. However, the critical current increases lin- 
early as the applied field becomes more negative. The 
maximum critical current [see Eq. (18)] is reached at the 
field [see Eq. (17)] when the critical current for vortex 
nucleation at a sharp inner corner becomes equal to the 
critical current for antivortex nucleation along one of the 
straight sides far from the outer corners. The behavior 
when the current makes right turns can be understood 
using simple symmetry arguments, as shown in Fig. 5(b). 

To confirm these theoretical predictions, we per- 
formed numerical simulations within the time-dependent 
Ginzburg-Landau (TDGL) theory. Figures 8, 12, and 
13 show results for the field dependence of the critical 
current in general agreement with the generic behavior 
predicted in Fig. 5. The predicted maximum enhance- 
ments of the critical current in negative applied fields 
calculated from Eq. (18) using the expressions for R and 
a given in Appendix A are Kcp/ Kcv(0) = 1.34, 1.53, and 
1.35, which are in reasonable agreement with the peak 
values of jc/jci^), in Figs. 8, 12, and 13 respectively. On 
the other hand, the fields at which the maximum en- 
hancements are predicted to occur, calculated from Eq. 
(17), Hmax/Hc2 = -0.012, -0.016, and - 0.014, are not 
in agreement with the positions of the peaks in Figs. 8, 
12, and 13, Hmax/Hc2 = -0.025, -0.03, and - 0.03, re- 
spectively. This discrepancy evidently arises from the 
quantitative inaccuracy of the London-model description 
and its handling of the vortex core on the length scale of 
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Another difference between the predictions of Sees. I- 

III and the TDGL simulations is that the peak of the crit- 
ical current vs field calculated using the London-model 
approach and shown in Fig. 5(a) resembles an inverted 
"V", whereas all the simulations show that this peak is 
rounded. The explanation of this difference is that the 
calculations of Sees. I-III were performed in the limit of 
very large W relative to ^, whereas the simulations of Sec. 

IV were carried out for W/^ equal to either 15.5 (Figs. 
8 and 12) or 13.5 (Fig. 13). Simulations for 90-degree 
turns with different values of W/^ (not shown here) have 
revealed that the peak in the critical current vs field be- 
comes sharper as the ratio W/S, increases. 

The above-described effects apply to all asymmetric ge- 
ometries of thin and narrow superconducting films. We 
thus predict that, upon application of a suitable per- 
pendicular magnetic field, the critical current in such 
patterns will depend strongly upon the direction of the 
current, allowing rectification of ac currents and produc- 
ing a diode effect. Another possible application of the 
field-induced enhancement of the critical current is in 
superconducting nanowire single-photon detectors. 
For the case of detectors with a two-dimensional mean- 
der layout (a "boustrophedonic" pattern), which have 
both left- and right-handed 180-degree turns, our results 
lead to the prediction that application of a perpendicular 
magnetic field would reduce the critical current. On the 
other hand, we predict that application of a negative ap- 
plied field would increase the critical current in detectors 
using thin and narrow superconducting strips in a spiral 
layout*^ where the current makes only left turns. Fab- 
rication of such spiral layouts, however, would be more 
difficult because of the need to electrically connect the 
center of the superconducting spiral through a hole in an 
insulating layer to a strip in a second superconducting 
layer. 
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Appendix A: Details 

1. Conformal mappings 

To calculate the functions needed for the calculation 
of the stream functions, sheet currents, and Gibbs free 



energies, we follow the approach given in Ref. 23 by first 
finding conformal mappings that map the upper half w- 
plane {w = u + iv) into the regions in the C-plane 
r + iy) corresponding to the films shown in Figs. 2-4. 



a. (.{w) for the sharp 90-degree turn 

The following mapping describes the sharp 90-degree 
turn shown in Fig. 2: 



dw 

aw) 



iW\/w — 1 



(Al) 



-KWs/w + 1 ' 

W 1 1 1 

— [tt -hi cosh Vw)-|-cosh V )],(A2) 

TT W 



and w(C), the inverse of C(^)i must be obtained numeri- 
cally. [This choice of w differs by 1 from that used in Eqs. 
(103) and (104) in Ref. 23.] Referring to Fig. 2, the letters 
A-F indicate mappings as follows: (A) x = y = +oo, 
w = — oo, (B) X = y = Q, w = —1, (C) x = +oo, y = 0, 
w = — e, (D) X = +O0, y = W, w = +e, (E) x = y = W, 
w = 1, and (F) x = W, y = +oo, w = +oo, where e is a 
positive infinitcssimal. 

Doing an expansion about the point E, we find that 
for ( a short distance 6 <^ W away from E along the 
diagonal between E and B [C = W{1 + i) - be"l'^\ the 
corresponding value of w is to lowest order 



w;(C) = l + i( 



37r5 \2/3 



(A3) 



h. C,{w) for the sharp rectangular 180-degree turnaround 

The following mapping describes the sharp rectangular 
180-degree turn shown in Fig. 3: 



dw 
C{w) 



(w2 - 1)^^2-/32' 



a 



(A4) 



(A5) 



where F{<l),k) and n((/), n, fc) are elliptic integrals of the 
first and third kind^^ of argument (/), modulus fc, and 
parameter n, and w{Q), the inverse of ^w), must be ob- 
tained numerically. Referring to Fig. 3, the letters A-J in- 
dicate mappings as follows: (A) x = W , y = 0, w = — oo, 
{B)x = W,y = W+h, w = -13, (C) x = -oo, y = W+h, 
w = — 1 — e, (D) X = — oo, y — h, w = (E) a; = 0, 

y = h, w = —a, (O) x = y = 0, w = 0, (F) a; = 0, 
y = —h, w = a, (G) x = — oo, y = —h, w = 1 — e, 
(H) X = -00, y = -W - h, w = 1 -\- e, (I) X = W, 
y = —W — h, w = (3, and (J) x = W , y = 0, w = +oo, 
where e is a positive infinitessimal. This mapping yields 



13 



three equations that determine the values of A, a, and /3 
for given values of W and h: 



W = 



W 



A 



(A6) 



- 1) 



-) 



)], (A7) 



h = A[K{-)-{l-a')U{a\-)], (A8) 

where K{k) and n(n, A;) are complete elliptic integrals 
of the first and third kind^^ of modulus k and parameter 
n. For example, solving Eqs. (A6)-(A8) for h/W = 0.15, 
as shown in Fig. 3, yields A = 0.609H^, a = 0.534, and 
/3 = 1.703. 

Doing an expansion about the point F for C along the 
diagonal between F and I, we find that for ^ a short 
distance S W away from F (C = -ih + ^e"'''/*) the 
corresponding value of w is to lowest order 

.(l-a2)(^2_ 2)1/3 2/3 

" = " + - 2aV3(^'-l)V3 (2^) • ^^'^ 



c- C(w) for the sharp 180-degree turnaround 

The following mapping describes the sharp 180-degree 
turnaround shown in Fig. 4: 



dC{w) _ .IW-JW 

dw IT 

Ciw) 
wiC) 



1 



W-i-Wtan-^ ( 



(w2 — l)-\/iir 



coth- 



2/ 7r(C-^)\ 
\ 2W J 



w 

/32csch2(- 



92)' 



(AlO) 
(All) 



where /? = cosh(7r/2), and the root in Eq. (A12) must 
be chosen such that Imu'(C) > 0. Referring to Fig. 4, the 
letters A-J indicate mappings as follows: (A) x = W, 
y = 0, w = -00, (B) X = W, y = W, w = -/?, (C) 
X = —00, y — W, w = —1 — e, (D) x — —00, y — e, 
w = — 1 + e, (O) X = y = Q, w = (G) x = —00, y = — e, 
w = 1 — e, (H) x = —00, y = —W, w = 1 + e, (I) x = W , 
y = —W, w = 13, and (J) x = W, y = 0, w = +00, where 
e is a positive infinitessimal. 

Doing an expansion about the origin O for ( along the 
X axis, we find that for ( a short distance S <^ W away 
from O (C = ^) the corresponding value of w is to lowest 
order 



coth 



2J 



1/2 



(A13) 



2. Stream function for the applied current 

To obtain the stream function for the applied current, 

we start with a complex potential Q^[w) describing the 
current flow in the tf-plane. The corresponding complex 
potential ^<^(C) = Qw{w{C)) describes the current flow 
in the ^-plane. For all values of x and y in the strip, 
Ki{x,y) = V X Si{x,y), where Si{x,y) = zSj^{x,y) is 
the stream function, and Siz{x,y) = ImC/j(x + iy), 
the imaginary part of GciO- In general, Kx{x,y) = 
dS{x,y)/dy and Ky{x,y) = —dS{x,y)/dx. 



a. Siz{x,y) for the sharp 90- degree turn 

For the sharp 90-degree turn, we start with the com- 
plex potential 



Qw{w) = Inw, 

TT 



(A14) 



which describes the flow of current I in the upper half 

w-plane from a source at w = 00 to a drain at w = 0. In 
the C-plane the same complex potential is 

GdO = GU^iO) = --H^iO], (A15) 

TT 

where w{C,) is the inverse of Ci^) given in Eq. (A2). The 
corresponding stream function Siz{x,y) is shown in Fig. 
2(a). 



b. Siz{x,y) for the sharp rectangular 180-degree turnaround 

For the sharp rectangular 180-degree turnaround, the 
complex potential 



^ . . I fw — l 
Qwiw) = - In I — 

TT V W + 1 



(A16) 



describes the flow of current / in the upper half w-plane 
from a source at w = 1 to a drain at w = — 1. In the 
C-plane the same complex potential is 

ec(C) = ^»MC)) = ^ln(^), (A17) 

where w{(^) is the inverse of Ci^^) given by Eq. (A5). The 
corresponding stream function Siz{x,y) is shown in Fig. 
3(a). 



c. Siz{x,y) for the sharp 180-degree turnaround 

For the sharp 180-degree turnaround, the complex po- 
tential is again given by Eq. (A16) or (A17), but with 
w{Q given in Eq. (A5). The corresponding stream func- 
tion Siz{x,y) is shown in Fig. 4(a). 



14 



3. Stream function for the field-induced current 



SHzix,y) for the sharp 90-degree turn 



Because V • Kh = 0, the sheet current induced by 
an applied field H ~ zH^ can be expressed as Kh = 

V X Sh, where Sh = zShz is the stream function. (This 
corresponds to the Meissner response of the sample.) 
This contribution to the total sheet-current density K 
must obey the London equation Kh — ~(2//ioA)[^ff + 
(0o/27r)V7ir], where Ah is the vector potential asso- 
ciated with the applied field {B = ^qH = V x Ah), 
and "fH is the corresponding phase of the order param- 
eter. There are many possible choices for Ah, but 
whatever choice is made, an expression for jh must 
be found that makes the gauge-invariant combination 
Ah + ((^o/27r)V7^f obey the same boundary conditions 
as does Kh, namely that (i) V • Kh = 0, (ii) Kh - n = 0, 
where h is an outward normal to the sample, and (iii) 
Kh carries no net current. 

We begin by separating Kh into two contributions, 
Kh = Ka + K^i, which obey Vxif^ = -{2/noK)B and 

V x K^ = 0, and we also separate the stream function 
into two contributions, Sh = Sa + S-^. With the choice 
A = yBz{x — W/2), we obtain 



KA = -y{2Hj^){x-W/2) 



(A18) 



and 



Sa = zSaz = z{H,/K)x{x - W). 



(A19) 



Although Ka has the properties required of Kh that (i) 
V • Ka = 0, (ii) Ka • n = along some of the boundaries 
shown in Figs. 2-4, and (iii) the integral of Ka across 
the strip is zero (the strip carries no net current) , it does 
not satisfy the requirements of zero normal component 
along all the boundaries of the sample. We therefore 
need to find K^ such that Kh = Ka + K-^ satisfies all 
the boundary conditions. 

We may represent the sheet-current density K^ = 
xK^x{x, y)+yKyy{x, y) by the analytic function K,~f{Q) = 
Kjx{x, y) — iK^y{x, y) of the complex variable C, = x+iy. 
The equations V x K^ = and V • K^ = are equiv- 
alent to the Cauchy-Riemann equations obeyed by ana- 
lytic functions. The corresponding complex potential is 
G-yiOi where A;^7(C) = dQ^{C,)/dC,, and the stream func- 
tion is S~f = zS-yz, where Sjz{x,y) — l'niQ^{x + iy). For 
each of the geometries shown in Figs. 2-4, the problem 
then reduces to finding the appropriate S-yz{x,y) that 
corrects for the failure of Saz to satisfy all the boundary 
conditions. 



For the sharp 90-degree turn, the stream function, 
Shz = Saz + Sjz, is given by 

SHzix,y) = ^{x{x-W) 

1, x(u)[x(u) -W] , ^ 

+ -Im / ^ ^ \ '-du}, (A20) 

n J_i w{x + iy) — u ) 

where x{u) = ReC(u -1- if), obtained from Eq. (A2)], is 

x{u) = (Pr/7r)[cos-i(-w) + cosh-i(l/|w|)], (A21) 

and w{x + iy) is obtained from the inverse of Eq. (A2). 
The resulting stream function is shown in Fig. 2(b) as a 
contour plot. 



b. SHz{x,y) for the sharp rectangular 180-degree 
turnaround 

For the sharp rectangular 180-degrec turnaround, the 
stream function, Shz = Saz + S^z, is given by 



SHz{x,y) = ^-^[x{x-W) 

<'^x{u) [x{u) — W\w{x+iy) 



A 

2 f^ ^W i^u) -wiw[x+iy) ^ 
TT Ja w'^{x+iy) — u'^ )' 



where x{u) = Rc(^(m + ie) is obtained from Eq. (A5) 
and w{x + iy) is obtained from the inverse of Eq. (A5). 
The resulting stream function is shown in Fig. 3(b) as a 
contour plot. 



c. SHz{x,y) for the sharp 180-degree turnaround 

For the sharp 180-degree turnaround, the stream func- 
tion, Shz = Saz + S^z, is given by 

SHz{x,y) = ^[x{x-W) 

''^x{u)[x{u) — W]w{x+iy) 
/o 



+ 



h4 

71" Jo 



2, ^. . 2 rf«|,(A23) 
w^yx+iy) — J 



where x(u) = ReC(ii -I- ie) is obtained from Eq. (All)] 
and w{x + iy) is obtained from Eq. (A12). The resulting 
stream function is shown in Fig. 4(b) as a contour plot. 



4. Applied sheet-current density near the sharp 
inner corners 

As suggested by the current crowding shown in Figs. 
2(a)-4(a), the sheet-current density generated by the ap- 
plied current diverges upon approaching the sharp in- 
ner corners. In the following we calculate Kig, the 



15 



counterclockwise component of the sheet-current den- 
sity at a distance 5 ((5 <C W) from the corner. We 
use the complex potentials given in Sec. A 2 to evaluate 
/Cc = dg/dC = K^- iKy. 



a. Kji for the sharp 90- degree turn 

Starting from Eqs. (A15) and (Al), and doing an ex- 
pansion about the point E using Eq. (A3), we find for ( 
a short distance S <^W away from E along the diagonal 
between E and B [( = W{1 + i) - Se^""/*] that the coun- 
terclockwise component of the sheet current density is to 
lowest order 



K 



IS 



K, 



/ 414^ \ 1/3 



(A24) 



V37r(57 

Integrating the Lorentz force (j)QKis, we find that the 
work done by the source of the current Ki in moving the 
vortex from the corner to the position 5 is 



Wi{5) = cf>aKi(^ 



1/3, 



(3(5) 



2/3 



(A25) 



b. Kis for the sharp rectangular 180-degree turnaround 

Starting from Eqs. (A17) and (A4), and doing an ex- 
pansion about the point F using Eq. (A9), we find for C, a 
short distance 5 -^W from F along the diagonal between 
F and I (C = —ih + ^e~"/^) that the counterclockwise 
component of the sheet current density is to lowest order 



v2U/3 



(^2-1)1/3 VSTraJ 



^ 2W y/3 



(A26) 



Integrating the Lorentz force (j)oKis, we find that the 
work done by the source of the current Kj in moving the 
vortex from the corner to the position 5 is 



Wi{6) = cPoK, 



(/32-a2)l/3 ^ 1/3 



(/32- 1)1/3 VTra 



)'( 



3(5 \ 2/3 

yJ ■ 



(A27) 



5. Field-induced sheet-current density near the 
sharp inner corners 

As suggested by the current crowding shown in Figs. 
2(b)-4(b), the sheet-current density induced by the ap- 
plied magnetic field (the Meissner current) diverges upon 
approaching the sharp inner corners. In the following we 
calculate Khs, the counterclockwise component of the 
sheet-current density at a distance 6 {S W) from the 
corner. We obtain Khs by evaluating = V x Sh, 
where Sh = zShz is one of the stream functions given 
in Sec. A 3. For each of the three cases, the divergence of 
Khs arises from the integrals given in Eqs. (A20), (A22), 
and (A23). 



a. Kh6 for the sharp 90-degree turn 

Starting from Eqs. (A20) and (Al), and doing an ex- 
pansion about the point E using Eq. (A3), we find for ( 
a short distance 6 W away from E along the diago- 
nal between E and B [C = W{1 + i) - de^''^^] that the 
counterclockwise component of the field-induced sheet- 
current density is to lowest order 



Khs = cr 



A KSttSJ 
where, with x{u) given by Eq. (A21), 



"=l^/-r 



x{u)[x{u)-W] 

(1 - U)2 



du 



8G 



(A30) 



0.742 (A31) 



and G = 0.915965... is Catalan's constant. Integrating 
the Lorentz force (PqKhs^ we find that the work done by 
the source of the applied field in moving the vortex 
from the corner to the position 5 is 



Wh{5) 



Poa- 



A 



(A32) 



c. Kis for the sharp 180-degree turnaround 



b. Kh6 for the sharp rectangular 180-degree turnaround 



Starting from Eqs. (A17) and (AlO), and doing an ex- 
pansion about the origin O using Eq. (A13), we find for C 
a short distance 5 -^W from O along the x axis {Q = 5) 
that the counterclockwise component of the sheet current 
density is to lowest order 



Kis = Ki{ — coi\i- 
\Tro 2 



1/2 



(A28) 



Integrating the Lorentz force (poKjg, wc find that the 
work done by the source of the current Kj in moving the 
vortex from the corner to the position i5 is 



/ W6 7r\ 
Wi{S) = 2ct>oKi(^—coth-j 



1/2 



(A29) 



Starting from Eqs. (A22) and (A4), and doing an ex- 
pansion about the point F using Eq. (A9), we find for ( a 
short distance (5 (/i ^ (5 <C W) from F along the diagonal 
between F and I (C = —ih -\- Se^^""^^) that the counter- 
clockwise component of the field-induced sheet-current 
density is to lowest order 



Khs - (T—f^ _ 1)1/3 I ' 



(A33) 



where, with x{u) obtained from Eq. (A5), 

(l-c^2) (a^ +u'^)x{u)[x{u)-W] 



a = 



(a2 - v?y 



du. (A34) 
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FIG. 14; Parameters a, /3, and the slope parameter o [Eq. 
(A34)] as functions of h/W for the sharp rectangular 180- 
degree turnaround shown in Fig. 3. The strip width is W , 
and the gap opening is A = 2h. 



6. Critical current for the nucleation of vortices 

We now use a London-model description of a nucleat- 
ing vortex to calculate the critical current for the nucle- 
ation of vortices. As in Ref. 23, the critical current is 
defined as the current that reduces the Gibbs free energy 

barrier at the nucleation site to zero, where the Gibbs 
free energy is the vortex self-energy Eseif{6), accounting 
for its interaction with all images, less the work Wi(^) 
done by the Lorentz forces Wl{S) in moving the vortex 
to the position S away from its nucleation site. In the 
presence of the applied field , this work term has con- 
tributions from both the applied current and the applied 
field: Wl{6) = Wi{6) + Wh{S). 



a. Critical current for the sharp 90-degree turn 



For example, solving Eqs. (A6)-(A8) for h/W = 0.15, 
as shown in Fig. 3, yields a = 0.534, /3 = 1.703, and 
a = 0.715. Similarly for h/W = 5/27 0.185, as shown 
in Fig. 13, a = 0.584, and 13 = 1.606, and a 0.721. As 
shown in Fig. 14, a varies slowly from 0.671 at h/W = 
[see Eq. (A37)] to 0.742 for large values of h/W > 1 [see 
Eq. (A31)]. 

Integrating the Lorentz force (j)oKHs, we find that the 
work done by the source of the applied field in moving 
the vortex from the corner to the position 6 is 



WHi6) 



, if^W^(^2_Q,2)l/3 1^ 1/3.3^.2/3 

A - 1)1/3 lyj • 



(A35) 



c. Kh6 for the sharp 180-degree turnaround 

Starting from Eqs. (A23) and (AlO) , and doing an ex- 
pansion about the origin O using Eq. (A13), we find for C 
a short distance S <^W from O along the x axis {( = d) 
that the counterclockwise component of the field-induced 
sheet-current density is to lowest order 



/^H. = a— (-coth-j 
where, with x{u) obtained from Eq. (All) 



1/2 



^ x{u)[x{u) - W] 



du = 0.671. 



(A36) 



(A37) 



Integrating the Lorentz force (poKjs, we find that the 

work done by the source of the current Kj in moving the 
vortex from the corner to the position 6 is 



iy.(^) = 2^o.'^(^coth|) 



1/2 



(ASS) 



For the sharp 90-degree turn, the Gibbs free energy 

becomes 



G 



In 



(f)-Wi{6)-WH{S), (A39) 



27T/ioA V e 



where the first term on the right-hand side was derived 
in Ref. 23 and the second and third terms are from Eqs. 
(A25) and (A32). Following the steps that led to Eq. (5), 
we obtain with Kj = K^v, 



Sc 



2TrnoMKi + a{WH^/A)] \4WJ 
6^/2^/3 = 1.49^, 



/ 37r \i/3 
I 1 (A40) 



K. 



where .^'^(0) = KcsR, K^s is given by Eq. (5), 

3/<Ni/3 



(A41) 
(A42) 



(A43) 



and a = 0.742, obtained from Eq. (A31). See Fig. 5(a). 
For the dimensions of Fig. 7,W = 15.5$, R = 0.555, and 
a = 0.742. 



6. Critical current for the sharp rectangular 180-degree 
turnaround 



For the sharp rectangular 180-degree turnaround, the 
Gibbs free energy becomes 

^=2|,A^^(f)-^^(^^-^-^^)' 

where the first term on the right-hand side was derived 
in Ref. 23 and the second and third terms are from Eqs. 
(A27) and (A35). Following the steps as in Sec. A 6 a, we 
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obtain Kcv{Hz) as in Eq. (A42), where -K'ct(O) = KcsR, 
Kcs is given by Eq. (5), 

^~ 2(/32-a2)i/3UvKJ ^^^^> 

and a must be obtained numerically from Eq. (A31) for 
a given ratio h/W. For the dimensions of Fig. 13, W = 
13.5^, h = 2.5^, a = 0.584, ^ = 1.606, i? = 0.545, and 
cr = 0.721. 



c. Critical current for the sharp 180-degree turnaround 

For the sharp 180-degree turnaround, the Gibbs free 

energy becomes 



where the first term on the right-hand side was derived 
in Ref. 23 and the second and third terms are from Eqs. 
(A29) and (A38). Following the steps that led to Eq. (5), 
we obtain with Ki = Kcv, 

1/2 ^ ^o[7rtanh(7r/2)/iy]V^ 
^ 27rfioA[Ki + a{WH,/A)]' ^ ' 

5^ = = 1.85^, (A48) 

and we obtain Kcv (Hz) as in Eq. (A42), where Kcv{0) = 
KcsR, is given by Eq. (5), 

i?=(^tanh-) , (A49) 

and a = 0.671, obtained from Eq. (A37). For the dimen- 
sions of Fig. 12, which we treat as a sharp 180-degree 
turnaround, W = 15.5^, R = 0.431, and a = 0.671. 
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